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Abstract 



The q— difference analog of the classical ladder operators is derived 
for those orthogonal polynomials arising from a class of indeterminate 



►^ , moments problem. 

m 

tJ- ■ 1 Introduction 

^. . 

This work is a follow up to our work [3] where we derived raising and lowering 
. operators for polynomials orthogonal with respect to absolutely continuous 

f"^ I measures /u under certain smoothness assumptions of //'. This approach 

goes back to [2], [3], and [TI]. The raising and lowering operators derived 
in these references are differential operators. It was later realized that a 
similar theory exists for polynomials orthogonal with respect to a measure 
C^ ' with masses at the union of at most two geometric progressions, {ag", 6g"}, 

for some q £ (0, 1), ^. The corresponding theory for difference operators is 
in [9]. This material is reproduced in [8]. The raising and lowering operators 
involve two functions An{x) and Bn{x) which satisfy certain recurrence re- 
lations. In the case of differential operators we have demonstrated that the 
knowledge of An{x) and Bn{x) determines the polynomials uniquely in the 
cases of Hermite, Laguerre, and Jacobi polynomials, see [5]. This is done 
through recovering the properties of the polynomials including the three 
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term recurrence relation which generates the polynomials. This work shows 
that the corresponding functions determine the polynomials in the cases of 
Stieltjes-Wigert and g-Laguerre polynomials. 

The orthogonal polynomials which arise from indeterminate moment 
problems have discrete and absolutely continuous orthogonality measures 
[1]. In many instances it is more convenient to work with absolutely contin- 
uous measures [8l, Chapter 21]. 

In this work we derive raising and lowering operators for polynomials 
orthogonal with respect to absolutely continuous measures. We shall assume 
that {Pn{x)} are monic orthogonal polynomials, so that 

(1.1) / w{x)Pjn{x)Pn{x)dx = CnSm,n- 

Jo 
A weight function w leads to a potential u defined by 



(1.2) u{x) 



Dq-iw{x) 



w{x) 
where Dq is the g-difference operator 

f{x) - f{qx) 



(1-3) iDJ)ix) 



X — qx 



Every monic sequence of orthogonal polynomials satisfies a three term re- 
currence relation of the form 

(1.4) (X - an)Pnix) = Pn+lix) + /3„P„_l(x). 

We also write the monic polynomials Pn{x) as follows: 

P„(x)=x" + pi(n)x"-i + ... 

and it follows immediately from the three term recurrence relations (1.4) 
that 

(1.5) an = Pi{n)-pi{n+l). 

A main result of this work is the following theorem. 
Theorem 1.1. Let 

(1.6) Ar,{x) := 1 r ^M-^(^) p^(^)p^(^/g) ^(y)rfy^ 

Cn Jo qx-y 

(1.7) B^{x) := -^ r ^(g^)-^(^) p„(^)p„„^(^/g)^(^)rf^. 

Cn-i Jo Qx-y 



Then we have the lowering relation 

(1.8) DgPnix) = (3nAn{x)Pn-i{x) - S„(x)P„(x). 

Theorem 1.1 will be proved in §2 along with the difference equations 
satisfied by A„(x) and i?„(x); 

n 

(1.9) Bn+i{x) + Bn{x) = {x- an)An{x) + x{q - l)^^j(x) -u{qx), 

j=0 

(1.10) 

1 + (X - an)Bn+l{x) - {qx - an)Bn{x) = (3n+lAn+i{x) - /?„^„_i(x). 

The identities (|1.9p - (|l.lU|) will be referred to as the supplementary condi- 
tions. 

Theorem 1.1 is the (/-analogue of 

P;l(x) = (3nAn{x)Pn^l{x) - 5„(x)P„(x) 

of [1] . See also [5j for a derivation of the supplementary conditions for the 
q = 1 case. 
Let 

(1.11) Lr,n:=Dg + Bn{x), 
be the lowering operator. Thus (jl.Sp is 

(1.12) Li^nPnix) = PnAn{x)Pn-l{x). 

The raising operator can be found as follows: First replace /?„P„_i(x) in 
(1.8) by (x — an)Pn{x) — Pn+i{x), using (1.4), and then replace 

-(x - an)An{x) + Bn{x) 

by 

n 
x{q- l)^^j(x) -u{qx) - P„+i(x), 
j=0 

using (1.9). An easy computation shows that, 

DqPn{x) = Bn+l{x)+u{qx)-x{q-l)^Aj{x) Pn{x)-An{x)Pn+l{x). 



With the replacement of n by n — 1 in above equation, the raising operator 
is, 

n-l 

L2,n ■= Dq + x{q - 1) X] Aj{x) - u{qx) - Bn{x), 
j=0 

and 

L2,nPn~l{x) = -An-i{x)Pn{x). 

It is useful to recall the following analogue of the product rule 

(1.13) Dq{f{x)g{x)) = {DJ{x))g{x) + f{xq)Dqg{x). 

The following lemma, whose proof is a calculus exercise, will be used in the 
proofs of our main results. 

Lemma 1.2. // the integrals 

dx f°° dx 

f{x)g{x) — , / f{x)g{qx) — , 
X Jo x 

exist then the following q-analogue of integration by parts holds 

(1.14) / f{x)Dqg{x)dx = — g{x)Dq-.f{x)dx. 
Jo Q Jo 

An immediate consequence of Lemma 1.2 and (jl.ip is 

(1.15) / u{y)Pn{y)Pn{y/q)wiy)dy = 0. 
Jo 

We also have 

f°° 1 - g"+^ 

(1.16) / u{y)Pn+i{y)Pn{y/q)w{y)dy = — gCn, 

Jo ^ — (1 

which follows from (fTTSl) . (fOll . (fLTil) . and the fact that 

1 — rt" 

Z?„x" = ^x"-^ 

l-q 



2 Proofs 

We shall need the fact [l3], [8] 

(2.1) Cn = Col3l(32 ■ . . Pn, 

and the ChristofFel-Darboux identity [13], [S] 



(2.2) 



n-l 

E 

fc=0 



Pfc(x)Pfc(y) P„(x)P„_i(y) - P„(2/)P„_i(x) 



a 



Cn-i(2;-y) 



Proof of Theorem 1.1. Let DqPn{x) = X]fc=o '^n,kPkix). Then 

Cfcc„,fc = / Pk{y)w{y)DqPn{y)dy. 
Jo 



Applying Lemma 1.2, p.l3p . we see that 

/■oo 

gCfcCn,fc = -/ P„(2/) [(i:'q-iPfc(y))u;(y) + Pfc(y/g)i:'g-iw(y)] rfy 



JO 

oo 



Pn{y)Pk{y/q) 



Dq-^w{y) 
w{y) 



w{y)dy 



where the orthogonality property was used in the last step. The definition 
of u dLlD yields 



Pn{y)Pk{y/q)u{y)w{y)dy 

Pn{y)Pk{y / q){u{qx) -u{y))w{y)dy, 



qCkCn,k 



where we again used the orthogonality property in the last step. Therefore 
by the Christoffel-Darboux identity (j2.2p 



1 



DgPnix) 



Cn-i Jo '" ■ qx-y 

and (12. ip . the theorem follows. 

Proo/ o/ ([L9]). It is clear that 

Bn+i{x) + Bn{x) 

u{qx) - u{y) 



D 



Qn{qx - y) 



[Pn+l{y)Pn{y/q) + PnPn{y)Pn-l{y/q)\ w{y)dy 



h + I'. 



2i 



where 

Qn Jo qx-y 

h ■■= J_ r ^(g^) - ^(y) [P^^,(y)p^(y/q) - P4y)P^_,,{y/q)] w{y)dy, 

Cn Jo qx-y 

after finPn^i{y/q) is replaced by {y/q - an)Pn{y/q) - Pn+iiy/q)- It is easy 
to see that Ji is given by 

h = {x- an)An{x) - —— / {u{qx) - u{y))Pn{y)Pn{y/q)w{y)dy 

^nQ Jo 

= {x- an)Anix) - q'^^'^uiqx), 
where (ll.lSj) and the fact that 

(2.3) Pj{y/q) = q^^ Pj{y) + lower degree terms 

were used. To evaluate I2 first note that ()2.3p implies 



r-oo 

(2.4) / P,{y)Pj{y/q)w{y)dy = Q^q-K 

Jo 

Next we apply the Christoffel-Darboux formula to 

Pn+i{y)Pn{y/q) - Pn{y)Pn+i{y/q), 

and replace y — y/q by {y — qx + qx){l — l/q). Therefore we see that 

h = x{q - 1) X; M^) +^-^^ r Hq^) - u{y)] J2 ^^^^^^^^^^w{y)dy 

=x(,-i)t^.(^)+^n(,x)t.---+^ri: '^^'^'^?^'^'^ ^.---(^)^^- 

,=0 1 j=0 ^ -^0 j.^0 ^^ 

Thus 

1-q , ,l-g~"-i 



h = x{q-l)S2Aj{x)^ u{qx)—- 



j=0 ^ 



Simplifying /i + I2 we establish ()1.9p . D 



Proof of (jl.lOp . From the definition of Bn{x) we see tliat 

[x - an)Bn+i{x) - {qx - an)Bn{x) 
' , Mqx) -u{y) 
qx-y 

"" " "-) Pn+l{y)Pn{y/q) - (^^^) Pn{y)Pn-liy/q) 



Cn 

w{y) [u{qx) - u{y)] 



dy 



Cn V? qx-y ' 



1 ^i + ^^ii)p^(y)p^_^(y/^)] 



dy 



C„-i V QX-y 

= --^ / w{y)u{y)Pn+i{y)Pn{y/q)dy 
CnQ Jo 

+1 r wiy f^"^""^ ~ ^^^^ (y/g - an)Pn{y/q)Pn+i{y)dy 
Cn Jo qx-y 

+7 / w{y)u{y)Pn{y)Pn-i{y/q)dy 

Sn-l JO 

1 r^(y)!^M^^p^(y)p„_^(y/^)(y_«„)rfy 



Cn-1 Jo qx-y 

1 /■°° 
= --— w{y)u{y)Pn+i{y)Pn{y/q)dy 
Cnq Jo 

Cn Jo ga; - y 

1 f°° 
+7 / w{y)u{y)Pn{y)Pn-i{y/q)dy 

^n-l Jo 

1 r^y)!^M^^(p„^,(y)+;3„p„_^(y))p„_,(y/^)dy. 



Cn-1 Jo qx-y 

The result follows after some simplifications using (|1.16p . D 

3 Stieltjest-Wigert polynomials 

The computations in this and the next section will show that An{x) and 
Bn{x) are rational in x and consequently (1.9) and (1.10), the supplemen- 
tary conditions which are identities valid for all x would be particularly 
useful for the recovery of the recurrence coefficients. As we shall see, sys- 
tems of apparently non-linear difference equations generated by equating 



the residues on both sides of (1.9) and (1.10), can be solved explicitly which 
ultimately determine the recurrence coefficients. 

This is example of an indeterminate moment problem associated with 
the log-normal density. See [6] for a discussion of the associated moment 
problem. We take the weight to be 

w{x) = cex.p[{lnx)'^/{2lnq)], 0<x<oo, < q < I, 

where c is a positive constant which will not appear in subsequent compu- 
tations. 

An easy calculation shows that 



q 


(^ Vq 


l-q 

Rn 

X2' 


\x x'^ 

1 r 


Cn(l ■ 

rn 

X2 


- q)V^ Jo 
l-q"! 
1 — q X 
1 



u{x) 

X- 

1 f°° dy 

where Rn := ——. :-^ / Pn{y)Pn{y/q)w{y) — , 

Cn(l - q)y/q Jo y 

r-_ 1 _ „" 1 

Bnix) 

1 /■°° dy 

where r„ := — / Pn{y)Pn-i{y / q)w{y) — . 

C,n-l^/q[l -q) Jo y 

From the supplementary conditions, p.9|) and p.lU|) 

n+l _|_ „n _ 2 1 

(3.1) ^__±y = i?„ + (g_i)5„ 



l-q 1 

1 



(3.2) r„+i + r„ = -a„i?„ 



^(1 - q) 

l_qn+l i_qn 

(3.3) = On—; an- \- Vn+l - qTn 

l-q l-q 

(3.4) a„(r„-r„+i) = (3n+lRn+l - PnRn-l, 

where Sn '■= Yll=o -^i' ^^^ 

1 !rMy)/y dy _ i 



(3.5) i?i 







^/qiX-q) J^w{y)dy 1 



A difference equation satisfied by Rn is found by subtracting (3.1) at "n — 1" 
from the same at " n" ; 



(3.6) qRn-Rn-i = -{l + q)q 



n-l 



and since the "integrating factor" is q ", the unique solution is 

n 

(3.7) R, ^ 



n — . 

1-q 

Note that (3.3) simphfies to 

(3.8) - anQ^^ = r„+i - qr^. 

Multiplying (3.2) by 1 — g', together with i?n(l — q) = q^ and (3.8) one finds 

(3.9) r„-gr„+i = — , 

v 9 

and since the "integrating factor" for this is g", the unique solution subject 
to ro = 0, is 

1 -a"" 

(3.10) Tn 



which with (3.8) immediately gives, 

-n 

(3.11) a„ = i^ (g-«-i + g-« - 1) . 

Multiply (3.4) by Rn and replace anRn with (3.2), we find the resulting first 
order difference equation 

' T) -I— 1 10 ' Tt \ 

r^ -ZTTZ r 1 = (3n+lRn+lRn — PnRnRn-1, 



where the solution with the initial conditions tq = /3o = is 

(3.12) rl - ?" = PnRnRn-l. 

\/9(i - q) 

This expresses /?„ in terms of the subsidiary quantities r„ and i?„. 



RnRn-l V \/^(l ~ ^ 

(3.13) = g-^"-g-3". 

In the next section we take a route for the computations of the recurrence 
coefficients which does not involve the determination of the analogous r„ and 

Rn- 



4 g-Laguerre polynomials 

This is also associated with an indeterminate moment problem at a level 
"higher" than the Stieltjest-Wigert polynomials, in the sense that when 
an appropriate limit of a parameter is taken, the g— Laguerre polynomials 
become the Stieltjest-Wigert polynomials. See pTO]. We refer the reader to 
[T2j for further information. 
We take the weight to be, 

(4.1) w(x) = ^, < X < oo, a > -1, < q <l, 

{-x;q)oo 

where 

oo 

{z;q)oo:='[lil-zq'') 
fc=0 

This weight leads to 

u{x] = h 



1 — q \ X X -\- q ^ 

u{qx) — u{y) 1 /?^" — 1 q~" 



qx — y l — q\ xy (x + l)(y + 

By definition, 

An{x) = -£l^ r p^{y)p,^{y/q)'!^dy 

Cn(l - q)x Jo y 



-a 



Pn{v)Pn{v/q)^^dy 



(4.2) 



C„(i-g)(x + i) 7o y + 1 



X {l-q){x + l) 

-a _-^ /-oo r \ 

Bn{x) = -^-— / Pn{y)Pn-l{y/q)^^dy 

(1 - q)C.n-ix Jo y 



P^{y)P^_,{y/q)^^dy 



Cn-iii-q){x + i) Jo y + q 



(4.3) 



rn g" ^Pi(ra) 



X X + 1 



where the evaluation of the second integrals in An{x) and Bn{x) will follow 
later. 



10 



Here is a computation of Rq. By definition, 

^ .^ g~"-l !^w{y)dy/y 

Q-" - 1 /(a) 1 



1-^ /(a + 1) 1-^ 
since 

I{a) := / —dy = — ■. . 

Jo {-y;Q)oo {q;q)oo smyra 

Also note the identity, 

1 1 1 



{-x;q)oo{x + q) q{l + x/q){-x;q)oo q{-x/q;q)oo 
Hence, 

°° w(v) f°° y" 

Pniy)Pniy/q)^dy = / Pn{y)Pn{y/q) . , , dy 

y + Q Jo q{-y/Q;Q)oo 

/>oo 

= 'Z" / Pn{qy)Pn{y)wiy)dy = g"+"Cn, 
Jo 

and the result for An{x) follows. Similarly, 

P^(y)P,,_,{y/q)^dy = r Pn{y)Pn-i{y/q) , f , dy 

y + Q Jo q{-y/Q;Q)oo 

/•oo 

= g" / Pn[qy)Pn-i{y)w{y)dy. 
Jo 

To complete the evaluation of the above integral, we note the identity, 

Pn{qy) = Pnto) + g"P„(y)-(z"P„(y) 

= q'^Pniy) + q^'y"" + (?"-ViNy""' + - - g" (y" + Pi(n)y"-^ + ...) 

= q^Pniy) + Pi(n)(g"-i - (?")P„-i(y) + ... 
Finally, 

oo 

Pn{qy)Pn^i{y)w{y)dy 



/•OD 

g° / {q'^Pniy) + Pi(n)(g"-i - q^Pn^iiy) + ...}P„-i(y)u;(y)dy 

JO 
--l')pi(n)g"+°Cn-i, 

11 



and the result for Bn{x) follows. 

It turns out that for the q—Laguerre weight, the supplementary condi- 
tions produce 6 difference equations in contrast with the 4 in the previous 
example. 

Now, by equating the residues for the simple poles at x = and x = —1 
in (|1.9p . we find, 



(4.4) r„+i + r„ = -OnRn z 



(4.5) pi(n + l)<7" + pi(n)g"-i = -l±^g" + L^!^ + ^, 

1-q 1-q 1-q 

respectively. We note here another identity involving i?„ only, by equating 
the constant terms of (|1.9p at x = oo. 



n i_ n+l 

(4.6) Rn - Y^- + {q- l)Sn - -YZ = 0' 

where Sn = ^21=0 ^j ■ ^ similar consideration on ()1.10p shows that 

(4.7) Onirn - rn+l) = Pn+lRn+l " PnRn-l 

(4.8) -(1 + a.)g>i(n + I) + {q + a„)g"-Vi(n) = ^"+'^ ^ _ ^"^ 

(4.9) r„+i - qrn = -g"a„ - 1. 

We use the fact that an = Pi{n) — Pi{n + 1) to rewrite (4.5) as a first 
order difference equation, 

Pi(n + l)g"+i - Pi(n)g"-i = g" + g"+i - (1 + q'") 

which has an "integrating factor" q'^~^ by inspection. Hence the above 
equation becomes, 

(4.10) Pi(n + l)g2" _ p^(n)g2"-2 = (i + g)g2n-i _ ^^ ^ g-")g"-\ 
and we find via a telescopic sum and the initial condition p^ (0) = 0, 

1 _ „2n+2 1 _ n+l 

pAn+l)q''' = ,/ , -(l+g-°)— ^. 
g(l - q) q{l - q) 

Therefore, 

(4.11) {l-q)piin) = -q+{l + q-)q-^+^-q-^^--+\ 

12 



and equation (jl.Sp gives 

(4.12) an = Pi(n) - pi(n + 1) = g-2n-i-"(l + q - q^+' - q^+»^'). 

At this stage Rn can be found from a difference equation obtained by sub- 
tracting (4.6) at "n" from the same at "n + 1", 



(4.13) qRn+i -Rn = g"+' - g", 

with the initial condition Rq = 1/(1 — g). Having now determined a„ in 
(4.12), r„ can be found from (4.9) with the initial condition vq = 0. We pro- 
ceed to the determination of /3„. Multiply (4.8) by the "integrating factor" 
g" and by 1 — g, 

- (1 + a„)(7'"(l - q)pi{n + I) + {q + a„)g2-i(l _ q)p^{n) 

(4.14) = /?„+ig2"+i-/?„<72n-i. 

The l.h.s. of the above is simplified to, 

(1 - q)q'^ (pi(n) - pi(n + 1)) - a„(l - q) (pi(n + l)q^^ - Pi(n)g2-i) 

= (1 - q)q^''an [1 - (pi(n + 1) - Pi{n)/q)] . 

With (4.11), the term pi(n -1- 1) — pi{n)/q simplifies to 

1 _ g-2n-a-l^ 

and consequently (4.14) reduces to the first order difference equation, with 
the initial condition Pq = 0, 

(4.15) /3n+i<z'"+' - Puq^""-' = (1 - 9)g~'-"«n. 
Taking a telescopic sum, and noting that Yll=o '^j = ~Pi('^)) 

n-l 

(4.16) /3„g2n-i ^ (i_^)g-i-«^„^ = _(i_^)^-i-ap^(^). 

Finally, 

/3n = -g-2"-"(l-(z)Pi(n) 

= -g-^"~"(-g + (1 + g"")g^"" - g-2n-a+l^ 

(4.17) = g-^"-2"+i(l-g")(l-g"+"). 



13 



It is interesting to note that in the computations of a„ and /3„, not ah the 6 
equations are required. We have used only (4.5) and (4.8). However, for an 
exphcit expression of the g— Ladder operators and therefore the determina- 
tion of r„ and Rn we need the equations (4.6), (4.9) and a„ in (4.12). 

We end with the remark that in the case of the classical Laguerre polyno- 
mials, a„ = 2n + l + a, and /3„ = J21=o ^j ~ "-("- + «); and this is analogous 
to (4.12) and (4.16)-(4.17), however, with appropriate modifications in the 
g— case. 
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